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Abstract 



A renormalization group (RG) theory of Goldstone mode singularities in the 
0(n > 1) -symmetric <p 4 model is discussed. This perturbative RG theory is claimed 
to be asymptotically exact, as regards the long-wave limit of the correlation functions, 
where it predicts a purely Gaussian behavior of the transverse correlation function. 
However, we show that the results of this theory are incorrect, and the Gaussian 
behavior originates from a rough error in mathematics. Other relevant perturbative 
theories are discussed, as well. 

1 Introduction 

Recently, our theory of grouping of Feynman diagrams has been extended to include 
the region below the critical point The two-point correlation function in Fourier 
representation has been studied, which depends on the wave vector k. Based on our 
method of systematic consideration of all diagrams, we have shown in [2] that the behavior 
of the longitudinal (Gii(k)) and transverse (Gj_(k)) correlation functions in 0(n > 1)- 
symmetric </? 4 model below the critical temperature T c is non-Gaussian. It means that 
Gj_(k) ~ a k~ x± and Gii(k) ~ bk~ x W with exponents d/2 < X± < 2 and Ay = 2X±-dis the 
long-wave (at k — » 0) solution of our equations at the spatial dimensionality 2 < d < 4. 
These results coincide also with a non-perturbative renormalization group (RG) analysis 
provided in our paper. 

It turns out, however, that some perturbative RG theory |12l II 3j . which is claimed to 
be exact by its founders, contradicts our results - it predicts a purely Gaussian behavior 
with Aj_ = 2. We have put some effort to clarify this question and have found that the 
Gaussian behavior in this theory originates from a rough error in mathematics. 

2 The critical analysis 

We consider a y? 4 model with the Hamiltonian 
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where the order parameter <p(x) is an n-component vector with components y?i(x), de- 
pending on the coordinate x, T is the temperature, h is an external field. The ip 4 model 
exhibits a nontrivial behavior in close vicinity, as well as below the critical temperature 
T c , if the order parameter is an n-component vector with n > 1. The related long-wave 
divergence of the longitudinal and transverse correlation functions (in Fourier represen- 
tation) at T < T c has been studied in [3] based on the hydrodynamical (Gaussian) ap- 
proximation. Essentially the same problem has been studied before in jjj in terms of 
the Gaussian spin-wave theory Later perturbative renormalization group (RG) stud- 
ies 0IZlEliniIiniiniII3IlllIIl support the Gaussian approximation. The RG method 
is claimed to be asymptotically exact. Hence, we disprove this statement. Our analysis 
in |2j predicts a non-Gaussian behavior, and we show by general physical arguments that 
it must be the true behavior to coincide with the know rigorous results for the classical 
XY model. 

According to the conventional believ [HI El E3 EH EEU El EE! E] , the transverse cor- 
relation function Gj_(k) diverges like k~ x± with Aj_ = 2 at k — > below T c for the systems 
with 0(n > 2) rotational symmetry. It corresponds to the Gy(k) ~ k d ~ 4 divergence of 
the longitudinal correlation function. Besides, the singular structure of the correlation 
functions is represented by an expansion in powers of k 4 ~ d and k d ~ 2 HUEj. Formally, 
our analysis in [2] is consistent with these results at Aj_ = 2, although in reality Aj_ < 2 
holds. Below we will show that Aj_ < 2 holds near two dimensions at n = 2. 

As usually accepted in lattice models, here we define that all the parameters of the 
normalized Hamiltonian H/T (pQ) are proportional to the inverse temperature 1/T. In this 
case ro is negative. The assumption that Gj_(k) ~ a(T) k~ 2 holds (with some temperature- 
dependent amplitude a(T)) in the stable region below the critical point, i. e., at T < T c /C, 
where C is an arbitrarily large constant, leads to a conclusion that the critical temperature 
T c continuously tends to zero at d — > 2 (supposed d > 2). Really, at Aj_ = 2 we have 

^ 2 ( x )) = ^E G *( k )^ M2 ( 2 ) 
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where Gj(k) = (| v?j(k) | 2 ) is the i-th component (one of which is longitudinal, other ones 
- transverse) of the correlation function, M is the magnetization, A is the upper cutoff 
of the wave vector magnitude, and S(d) is the area of unit sphere in d dimensions. Here 
G||(k) = Cr'|(k) + 5k : oM 2 V is the Fourier transform of the real-space correlation function 

((/?H (0)<^|| (x)), whereas Gy(k) denotes the Fourier transform of (y?||(0)^||(x)) — M 2 , where 
<^||(x) is the longitudinal component of the order-parameter field. Since the amplitude 
of the transverse fluctuations never can vanish at a finite temperature, Eq. ((2j) implies 
that the average (i^ 2 (x)) diverges at T = T c /C when d — * 2, if T c remains finite. Thus, 
we obtain an unphysical result unless the critical temperature T c and, therefore, a(T c /C) 
tend to zero at d — > 2. 

On the other hand, it is a rigorously stated fact ^3E1 that the classical 2D XY model 
undergoes the Kosterlitz-Thouless phase transition at a finite temperature T KT . It means 
that a certain structural order without the spontaneous magnetization exists within the 
temperature region T < T KT . There is a general tendency of disordering with decreasing 
the spatial dimensionality d, and not vice versa. Thus, since the structural order exists at 
T < T KT and d = 2, some kind of order necessarily exists also at T < T KT and d > 2. Since 
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the classical XY model undergoes the disorder — > long-range order phase transition at 
d > 2, this obviously is the long-range order. Thus, the critical temperature at d = 2 + e 
is T c > T KT ^ for an infinitesimal and positive e, and it drops to zero by a jump at 
d = 2 — £, as consistent with the rigorous consideration in |15j . 

The classical XY model belongs to the same universality class as the actual ip 4 model 
at n = 2, which means that both models become fully equivalent after a renormalization 
(a suitable renormalization has been discussed in 2 ). Thus, T c does not vanish at d — > 2 
(for d > 2) also in the c/? 4 model. In such a way, the assumption Gj_(k) ~ a(T) k~ 2 
leads to a contradiction. In the stable region T < T c /C, the Gaussian approximation 
Gj_(k) ~ l/(2c/c 2 ) makes sense at finite not too small values of k. The above contradiction 
means that the Gaussian approximation with Aj_ = 2 cannot be extended to k — * in 
vicinity of d = 2. The contradiction is removed only if Aj_ < 2 holds at d — > 2 in the actual 
case of n = 2. 

It has been stated in ^21 EEU that the essentially Gaussian result Aj_ = 2 of the 
perturbative RG theory should be exact. However, the underlying method is not rigorous: 
it is assumed without proof that the renormalized Hamiltonian has the form of the Landau- 
Ginzburg- Wilson expansion with few terms included, as in the original Hamiltonian ijT)). 
Besides the "exact" rusult is obtained by cutting (at one loop order) a purely formal 
divergent perturbation series, where the expansion parameters are in no sense small. One 
claims |12 (I13[ IT4*] that the one-loop-diagram approximation provides asymptotically exact 
results at ro — > — oo (which corresponds to mo — > oo or m — > oo in ^JEU) an d, according 
to the provided there renormalization group arguments, also in the long-wavelength limit 
k — > 0. However, some of these "exact" results are rather unphysical. In particular, we 
find from Eq. (3.6) in ^3] and from the formula (tt 2 ) = — QA/uq given in the line just 
below that 

(vr 2 (x)) = (iV-l) J ^ (3) 

holds, where vr(x) is the transverse (N — l)-component field. Eq. © represents a sense- 
less result, since (7r 2 (x)} given by this equation diverges at d — > 2. It is clear that (7r 2 (x)) 
cannot diverge in reality, as it follows from the Hamiltonian density (2.1) in !13 a (Hamil- 
tonian (Q) in our paper): any field configuration with diverging vr 2 (x) provides a divergent 
7r(x)-dependent term 

~ ~ I Vtt(x) I 2 +|(vr 2 (x)) 2 (4) 

in the Hamiltonian density and, therefore, gives no essential contribution to the statistical 
averages. The result (J3J) corresponds to a poor approximation where the second term in Q 
is neglected. In a surprising way, based on Ward identities, authors of JB] and related 
papers have lost all the purely transverse diagrams, generated by the term ^(7r 2 (x)) 2 , 
and stated that this is the exact result at tq — > — oo. These diagrams, evidently, cannot 
vanish due to tq — > — oo, simply, because they are independent of r$. According to [T2Hr5] . 
the actual transverse term appears to be hiden in a shifted longitudinal field s, which is 
considered as an independent Gaussian variable (cf. Eqs. (3.5) and (3.6) in 12 j). Obviously, 
this is the fatal trivial error which leads to the above discussed unphysical result (|3"j). since 
s is not an independent field in reality - the determinant of the transformation Jacobian 
(from 7r, s to it, s) is omitted in the relevant functional integrals! Such an introduction 
of the shifted field to obtain a formally Gaussian Hamiltonian does not make any sense, 
because the transformation Jacobian then must be included and the resulting model all 
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the same is not Gaussian. Since © comes from 

(vr 2 (x)) = (iV-l)(2vr)-' i y' G ± (k)d d k, (5) 

the unphysical divergence of (tt 2 (x)) means that the predicted Gaussian form of the trans- 
verse correlation function Gj_(k) is incorrect. Another aspect is that the method used 
in ^3 gives Aj_ = 2 also at n = 2 in contradiction with our previous discussion 
concerning the known rigorous results for the XY model. 

Our consideration does not contradict the conventional statement (see ^1] and refer- 
ences therein) that the Gaussian spin-wave theory [3] becomes exact at low temperatures, 
but only in the sense that it holds for any given nonzero k at T — > 0, and in the limit 
lim lim in particular. However, the actual limit of interest is k — > or, equally, lim lim. 

Therefore, it is impossible to make any rigorous conclusion regarding the exponent Aj_ (or 
any related exponent) based on the fact that the Gaussian spin-wave theory becomes exact 
at T — ► 0. One has to prove that the limits can be exchanged! This problem persists also 
in the classical treatment of the many-particle systems in [I] which, in essence, is based 
on the Gaussian spin wave theory at T — > 0. This treatment, evidently, is not exact, since 
it breaks down at d — > 2 for the two-component (n = 2) vector model (where T c remains 
finite and we fix the temperature < T < T c ) just like we have discussed already - the 
average (7r 2 (x)} is given by the integral (6.8) in 4 which diverges in this case (supposed 
(2Tr)- 3 d 3 k -» {2^)- d d d k). 

A slightly different perturbative RG approach has been developed in ^U] to analyze 
the nonlinear a model. In this case the modulus of y(x) is fixed which automatically 
removes the divergence of (7r 2 (x)). A finite external field h has been introduced there 
to make an expansion. The correlation functions have the power-like singularities of 
interest only at h = +0, which means that in this case we have to consider the limit 
lim lim, i. e., the limit h — > must be taken first at a fixed nonzero k (p in formulae 

used in The results in JOJ are not rigorous since the expansions used there are 

purely formal, i. e., they break down in this limit. Besides, the renormalization of jlUj 
is no more than an uncontrolled approximation scheme: contrary to the assumptions 
in ^0], it should be clear that the exact renormalization is a fairly nontrivial problem 
which cannot be reduced to a finding of only two renormalization constants. If, e. g., 
we make a real-space renormalization of the Heisenberg model, say, with the scaling 
factor s = 2, then the statistically averaged block-spins of the Kadanoff transformation 
(composed of s d original spins) do not have a fixed modulus - simply, the original model 
does not include the constraint | ip(x) \= const for the block averages. It means that the 
transformation with any finite s yields a Hamiltonian of form different from the original 
one, i. e., the original Hamiltonian with merely renormalized coupling constant can never 
be the fixed-point Hamiltonian. This is the common problem in the perturbative RG 
theory that the renormalized Hamiltonian is assumed (or even "proved" like in ^U]) to 
have the same form as the original one. This, however, appears to be wrong in view of 
the exact renormalization. 

In absence of rigorous proofs, theoretical predictions can be supported by appropriate 
empirical data. However, the experimental measurements of susceptibility x depending on 
field h in isotropous ferromagnets like high-purity polycrystalline Ni are incompatible 
with the prediction \ ~ h^ d ~ A ^ 2 of the "exact" RG theory. An explanation has been given 
in |5j that the divergence of the susceptibility is not observed in real crystals due to the 
symmetry breaking perturbations causing easy axes. However, this argument is not valid 
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for polycrystalline materials, since they are practically isotropous in a length scale which 
exceeds remarkably the size of grains. Thus, the theory should work at h — > 0, where the 
relevant length measure, i. e. the correlation length, diverges. 

3 Conclusions 

The "exact" RG theory of Goldstone modes developed in ^1 and related papers is 
not exact, it is unphysical and mathematically erroneous. The essential claims about the 
Gaussian character of the 0(n > l)-symmetric if 4 model below T c are based in ^2] on the 
fact that the Hamiltonian can be formally written in an apparently Gaussian form. The 
author, however, forgot to include the determinant of the transformation Jacobian in the 
relevant functional integrals, according to which the resulting model all the same is not 
Gaussian. Surprisingly, the experts in the perturbative RG theory have not mentioned 
this evidently rough error in their theory till now. Other versions of the perturbative RG 
theory do not have any rigorous mathematical basis, as well. Besides, the prediction of the 
"exact" RG theory that the susceptibility below T c diverges as /i -1 / 2 in three dimensions 
evidently disagree with the experimental measurements in isotropous ferromagnets. 
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